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Abstract. We consider Schrodinger equations and Fokker-Planck equations in one 
dimension, and study the low-energy asymptotic behavior of the Green function using 
a new method. In this method, the coefficient of the expansion in powers of the wave 
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remainder term can be analyzed on the basis of an expression in terms of transmission 
and reflection coefficients. This method is applicable to a wide variety of potentials 
which may not necessarily be finite as a; — ^ ±cx). 
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1. Introduction 




(1.1) 



(1.2) 



The Fokker-Planck equation (1.2) describes the diffusion process in an external potential 
V{x), which is related to the function f{x) in (1.2) by 



or V{x) = —2 J f{x) dx. The correspondence between equations (1.1) and (1.2) is given 
by the relations 0(a;) = e~^(^)/^'0(a;) and 



Here we study the Green function for (1.1) or (1.2). Our analysis is based on the 
expression of the Green function in terms of reflection coefficients which was derived 
in [2], and the method of asymptotic expansion for the reflection coefficients presented 
in [3]. In this method, the high- and low-energy expansions can be treated on an equal 
footing. The high-energy expansion of the Green function was discussed in a previous 
paper [4]. In the present paper, we study the low-energy expansion. 

The low-energy asymptotic behavior of the Green function and related functions 
has been studied for many years, and various methods have been proposed [5-15]. (Most 
of these methods are limited to the cases where |Vs(a;)|(ia; < oo, although there are 
some specific methods for other types of Vs.) In this paper, we use a method different 
from any of these previous works, and derive simple results for the expansion in powers 
of the wave number k to arbitrary order. There have been previously similar attempts 
at systematically calculating the expansion of the Green function to arbitrary order 
[9], but the formulae derived in this paper are new, and applicable to a larger class of 
potentials (including the cases where Va,[x) is not finite as a; — )> ±oo). These results are 
not only of theoretical interest, but also useful for practical calculations. 

We assume that the Fokker-Planck potential V{x) is a real- valued function which 
is piecewise continuously differentiablej . We allow V{x) to be either -|-oo, — oo or finite 

X The Fokker-Planck equation (1.2) is well defined even when V{x) has a jump discontinuity and f{x) 
has a delta function [1], although the corresponding Vs{x) does not make sense in such a case. A 
delta function in /(x) is interpreted as jump conditions for (j){x) and (f>'{x), in much the same way 
that a delta fimction in Vs{x) is interpreted as a jump condition for 4''{x). The jump conditions 
require that '•^-'^(a;) and (p'{x) — 2f{x)(j){x) be continuous. For example, if f{x) — c5{x — xq), then 
= e^'^^(a;o — 0) and (I)'{xq+0) = ^'{xq — O)- If /(a;) has ajump, then Vs{x) has a delta function, 
and ■4''{x) has a jump. The Green function studied in this paper is meaningful for the Fokker-Planck 
equation even when Vs does not make sense. 




(1.3) 



Vs{x)^f'{x) + f{x). 



(1.4) 



Low-energy asymptotic expansion of the Green function 



3 



as X — >■ +00, and similarly for x — )> —00, and we require that V{x) behave steadily 
and smoothly at spatial infinity. Specifically, we assume that V{x), f{x), and f'{x) are 
monotone for sufficiently large at both x — )■ +00 and x — )■ —00. We do not consider 
the cases where V{x) shows oscillatory or other indeterminate behavior as a; — )> ±00. 
Note that Vs(±oo) are either finite or +00. The cases Vs(±oo) = —00 are not considered 
here, since such Vs does not correspond to an appropriate Fokker-Planck potential. 

Our aim in this paper is to derive an expansion of the Green function in powers 
of k for Imk > 0. For such an expansion to be possible, it is necessary that V{x) 
either converge sufficiently rapidly or diverge sufficiently rapidly as \x\ — > ±00. Let us 
introduce the following classes of real-valued functions: 

F^"-* ~ I 5' J (-'- + kr")l5'(^)l dx < 00 for any finite a j> , (1.5a) 

= I 5* J (1 + |a;|")|5((a;)| dx < 00 for any finite a| , (1.56) 

where n is a nonnegative integer. We derive the low-energy expansion under the 
condition that 

V-V^e F^-^ or e'^ e F^^ or e F^^ (1.6) 

with some n and some finite constant Vi. (Here V ~ Vi, e^^, and mean V{x) — Vi, 
6"^^^^ and e^*-"^-* as functions of x.) This is essentially a condition on the behavior of 
V{x) as X — >■ —00, and the three cases in (1.6) correspond to the cases V{—oo) = Vi, 
V{—oo) — +00, and V{—oo) — —00, respectively. Similarly, concerning the behavior of 
V{x) as X ^ +00, it is required that 

y-^seFW or e-^eFW or e (1.7) 

with some m and some finite constant V2. The three cases in (1.7) correspond to the 
cases V{-\-oo) = V2, ^(+00) = -l-oo, and V{-\-oo) — —00. Since there are three cases for 
X — > —00 and three cases for x — > -|-oo, there are nine cases in all. But it is sufficient to 
consider only the following six cases: 

(i) V{-oo) = Vi, V{+oo) = V2, (ii) V{-oo) = V,, V{+oo) = +00, 

(iii) V'(-oo) = Vi, V{+oo) = -00, (iv) V{-oo) = +00, ^(+00) = +00, 
(v) V{—oo) = +00, V{+oo) = —00, (vi) V{—oo) = —00, V{-\-oo) = —00. 

We derive an expansion of the Green function for each of the six cases. (The detailed 
conditions for the validity of the expansion to order are given by (5.29) and (6.3)). 

In section 2 we review the expression of the Green function derived in [2], and make 
a comment on its application to the Schrodinger equation. In section 3, we explain the 
method of [3] with explicit calculations. Using this, we derive the expressions for the 
expansion of the Green function in sections 4 and 5. The remainder term is studied 
in section 6. In section 7, we discuss the special case where Vs{x) tends to at both 
X — > -|-oo and x — > — 00. Examples of calculations are given in section 8. 
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2. Reflection coefficients and the Green function 



-Vsix) + e 



Gs{x,y;k)^5{x-y) (2.1) 



Let Gs{x, y; k) denote the Green function§ for the Schrodinger equation (1.1), satisfying 

dx"^ 

with the boundary conditions Gq,{x^ y;k)-^0 as 1/| — )■ oo for Im A; > 0. For Im k = 0, 
we define|| Gs{x, y; k) = lira^^o Gs{x, y;k + ie). Since G^{x, y; k) = Gs{y, x; k), without 
loss of generahty we assume x > y. 

In this paper, we use the expression of the Green function in terms of reflection 
coefficients for semi-infinite intervals. First, let us define the transmission and reflection 
coefficients for finite intervals. For arbitrary a and b {a <b), we define 

{V(a) (x < a) 
Vix) ia<x<b) f{x) = -\^V{x), (2.2) 
V{b) {b<x), 

and consider the Fokker-Planck equation (1.2) with f{x) replaced by f{x). Since 
f{x) — ior X < a and x > b, this equation has two solutions of the forms 

^[VQ,)-V[a)]/1^^-ik(x-a) < ^) 

0i(x) = < ., , ., , (2.3a) 

(The factor ^^'•''^1/^ in front of r is necessary in order to make the same coefficient 

r appear in both (2.3a) and (2.3b), since 0i and 02 are solutions of the Fokker-Planck 
equation. This factor disappears if we rewrite (2.3) in terms of ^pl{x) = e^^^^/^0i(x) 
and i^2{x) = e^(^)/^02(2;), which are solutions of the Schrodinger equation.) This 
defines the transmission coefficient r, the right reflection coefficient Rr, and the left 
reflection coefficient Ri for the interval (a, 6). We write them as T{b,a;k), Rr{b,a;k), 
and Ri{b,a]k). The reflection coefficients for semi-infinite intervals are defined as the 
limit a — )■ —oo of Rr{b,a;k) and the limit b — )■ -|-oo of Ri{b,a;k). When k is a real 
number, it may happen that these limits do not exist. In such cases, we define 

Rr{b,—oo;k)=lim lim Rr{b,a; k + ie), (2.4a) 

e4,0 a->— oo 

Ri{oo,a;k) = lim hm Ri{b,a;k-\-i€). (2-4&) 
In appendix G, it is shown that these limits indeed exist. 

§ The Green function for the Fokker-Planck equation (1.2) is GF{x,y; k) = e^[^(^'^^(^)l/^Gs(a;, y; k). 
When Vs does not make sense (see the footnote in the introduction), we need to define Gp first as the 
Green function for the Fokker-Planck equation, and then define Gs = e[^^^'~^^^^l/^GF- 
II This Gs(fc) becomes infinite where is an eigenvalue of the Schrodinger operator, and whc;re k 
corresponds to a half-bound state. Elsewhere, this definition of Gs{k) makes sense for real k, even on 
the continuous spectrum of the Schrodinger operator. 
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S^i^x^ /c) 



Si{x,k) 



Ri{oo; x; k) 
1 + Ri{oo, x; k) ' 



1 + Rr{x, — oo; k) ' 
S{x, k) = Sr{x, k) + Si{x, k). 
Then the Green function can be expressed in terms of this function S as [2] 



Gs{x,y;k) = 



1 



exp 



ik{x — y) — ik I S{z, k) dz 

J V 



(2.5) 
(2.6) 

(2.7) 



2ik^[l-S{xM[^-S{yM 
The function S{x, k) can also be expressed in terms of reflection coefficients for 
the Schrodinger equation. Let us consider the Schrodinger equation with the truncated 
potential which is set to be zero outside the interval (a, h): 





[X < a] 



dx' 



;ip{x) + Vs{x)i>{x) = k'^ip{x), 



Vs{x) ^ { Vs{x) 




{a<x<h) 
{b < x). 



(2.8) 



[x < a) 



This equation has two solutions of the forms 

The transmission coefficient T^{b,a;k) and the reflection coefficients Rf{b,a;k), 
Ri{b,a;k) are thus defined for the Schrodinger equation. We define the Rf and Rf 
for semi-infinite intervals in the same way as (2.4). It can be shown that the Sr and Si 
defined by (2.5) are expressed in terms of Rf and Rf as 

Rf{x,-oo) f{x) 



(2.10) 



l + Rf{x,-oo) 2ik' -'v-'-/ l + Rf{oo,x) 2ik ' 

(See appendix A for a proof.) The f{x) in (2.10) cancels out when we substitute (2.10) 
into (2.6), and so S takes the same form as the expression in terms of Rr and Rf. 

Rf{x,-oo) Rf {oo,x) 



S{x, k) 



(2.11) 



1 + R^{x,-oo) 1 + Rf {oo,x)' 

(Incidentally, note that the right-hand side of (2.11) can also be written in terms of the 
Weyl-Titchmarsh m- function.) 

3. Formulae for generalized reflection coefRcients 

Let us first make some definition. We define, for n = 1, 2, 3, . . . and — oo < a <b < +00, 

(3.1) 



[(7i, (72, ... , an]l =/■■■/ dzi---dznexp 

J J a<zi<Z2<—<z„<b 
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where each aj is either +1 or —1. The integrals mean dzi J^^ dz2 J^^ dz^ ■ ■ ■ ^ dz„ 
When Vi = hmj;^_oo V{x) is finite, we use the notation^ 

(-1, (72, ... , ar,t^ = [-1, (72, ... , a^t^ - [+1, a^, . . . , (7„]i^ 



2e 



-Vi 



0O<Zl<Z2<---<Zn<b 



dzi ■ ■ ■ dzn sinh [V^ — ^^(2:1)] exp 



Y,^,V{z,) 

■j=2 



. (3.2a) 



For n = 1, this means (— 1]^ 



2e ^i/^^sinh[yi — V{z)]dz. Similarly, when 



Vo = lim^^ 



V{x) is finite. 



[cri,...,cr„_i,+l]~ 



■n-l 

■j=l 



{3.2b) 



— 2e / • • • / dzi ■ ■ ■ dZn sinh [V2 — V{zn)] exp 

J J a<Zl<Z2< — <Zn<00 

The conditions for the existence of these integrals will be discussed later. 

In the formalism of [3], we deal with the generalized scattering coefficients, which 
are defined with an additional variable W as 

Rr{x,y;k)-ax, W) 



Rr(x,y;W;k) 



Ri{x,y;W;k) = Ri{x,y;k) + 
f{x,y]W]k) 



l-ax,W)Rrix,y; k)' 

i{x,W) T^{x,y;k) 



l-C{x,W)Rr{x,y; k)' 
-i{x,W) T{x,y;k) 



l-Cix,W)Rrix,y; k)' 



where 



^{x, W) = tanh 



W - V{x) 



-f{x, VT) = V 1 - = sech 



W - V{x) 



(3.3a) 
(3.36) 
(3.3c) 

(3.4) 



The original scattering coefficients Rr, Ri and r are recovered from i?^, Ri and r by 
setting W = V{x). We define the operator which acts on functions of x and VT, as 



Cg{x,W) = 2 



cosh[iy - V{z)\ + SYnh.\W - Vi^)]^ ) 9{z. dz. 



(3.5) 



The low-energy expansion of Rr for semi-infinite intervals was studied in [3]. 
According to the formulae derived there, we have 

Rr{x, -00; W] A;) = fo + ikfi + {ik)^f2 H h {ik)^fN + Pn, 



where"^ 



Tq — — tanh 



W -V{-oo) 



fn{x,W)^C^[fQ + i{x,W)] (n>l), 

f'^{x.z\W\k: 



PN = {ikf^' 



r'N+i{z, A) dz, 



,l-Rf{x,z;W;k) 

^ The relation with the notation used in [3] is e^^ 1, • • •]'Loo = (i; ' • •]'' 
+ For real fc, the integral in (c 
in the integrand, if necessary. 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



For real fc, the integral in (3.9) should be understood as lime4,o linij/-).-oo Jy with k replaced hy k + ie 
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The integer N > is arbitrary as long as fo,...,fN are finite. Let us derive the 
expressions for f„ and in terms of the integrals defined by (3.1) and (3.2). We consider 
the three cases, (a) V{—oo) = Vi (finite), (b) V{—oo) = +oo, and (c) V{—oo) = — oo. 
We let the superscripts a, b and c stand for the cases (a), (b), and (c), respectively 

(a) V(-oo) = Vi. 
In this case we have 

R,^r^ + ikfl + {ikffl + • • • + (ifc)^f^ + p%, (3.11) 

W -V^ 

f^ = -tanh — ^. (3.12) 

We can calculate by substituting (3.12) into (3.8). The details of the calculation are 
given in appendix B. As a result, we obtain, for n > 1, 

<(^, = E ^-i>-.>.--n-i (^) (-1' ^1' ^2, . . . , (3.13) 

{fl-'-CTn-l} 

where ^ i = V"^ _ , , Y]^ _ , i • • • V'^ _ , , , and 

cxD n— 1 

^.„...,.„_.(W^) ^ 2E(-l)"^+'^^S.,.„-.e^'""^^'^^e^'""^^'^. A ^ E'^^' (3-14) 

(3.15) 



m=l i=l 

n-1 



Substituting (3.13) into (3.9), the expression for the remainder term is obtained as 



/X —2 



{(Tl---(TJV} 



(3.16) 

(Here A = ^T^=\^i-) (3-16), f and stand for f{x^z\W\k^ and ^^(x, 1^; A;), 
respectively. (If = 0, the expression (— 1, . . . , (Tjv-i]1oo (3.16) is replaced 

by 2e-^isinh[Vi - V{z)\. Then A = and p(™) = 1.) For each A^, the right-hand 
side of (3.16) can be written in a compact form without the infinite sum over m. (See 
appendix B.) For example, defining C, = tanh{[Vi — K(2;)]/2}. we can write and pi as 

(3.17) 

Equation (3.11) makes sense if and only if fg, fj, . . . , all exist as finite quantities. 
(By construction, the remainder term is automatically finite if all are finite, since 
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Rr itself is finite.) Since f° (n > 1) has the form of (3.13), it is necessary that V{x) 
tend to Vi fast enough, for otherwise the integrals defined by (3.2) do not exist. We can 
show that is finite if V — Vi E (See appendix C.) Therefore, fg, fj, . . . , are 

all finite, and hence the expression (3.11) makes sense, if V — Vi E -^jv-i- 

(b) y(-oo) = +oo. 

Next, we consider the case V{—oo) — +oo. We write 

Rr = f'o + ^kf\ + {ikffl + ■■■ + iikff% + p%. (3.18) 

Unlike (3.12), the first term Tq, which is obtained from (3.7), is independent of W: 

= 1. (3.19) 

The higher-order coefficients are obtained by substituting (3.19) into (3.8). The 
calculation is easier in this case (see appendix B), and we obtain 

r'ni^. W)- Yl 2P«.....,..-.e^'~^^^ [-1' '^i' ^2, . . . , a^-it^, (3.20) 

{(Tl-"(Tn-l} 

/X —2 
dz [-1, <7i, <72, . . . , <7^r-i]loo e""^^'^ 

xP«...,..e(--)-(^)(i±|)"'. (3.21) 



(As before, A = '^i (3-20) and A = XlHi in (3-21).) 

For (3.20) to exist as a finite quantity, it is necessary that e~^^^^ tend to fast 
enough as x — ?> — oo. It can be shown that is finite if E F^Z}i (see appendix C). 
So the expression (3.18) makes sense if e -Fjv-i- 

(c) V{~oo) = — oo. 

The expressions for the case V{—oo) = —oo can be obtained in the same way. We have 
Rr = f'o + ikfl + {ikyf'2 + --- + {ikff''^ + p%, f(; = -l, (3.22) 

f^ix, W)^- J2 2Pi:!.„...,.„_,e-(^-^)^ [+1, -a,, -a,, (3.23) 
4. The low-energy expansion of Sr 

The function Sr defined by (2.5) can be extracted from Rr as 
S(x k)- ]im^-±^- hm ^i^^W) + Rr{x,-oo;W;k) 

br{x, k) - lym^ ^ _ - ^hm^ ^ _ ^^^^^ ^^^^ . (4.1j 

Since ^(x, W) ~ -l + 2e^-^(^) and 1 - [C(x, W)]'^ ~ ie^-n^) ^ q, we may write 

Sr(x, k)-- = \ lim e-^+^(^) [-1 + Rr(x, -oo; W; k)] . (4.2) 
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(We have moved the | to the left-hand side for convenience.) As in the last section, we 
consider the three cases (a), (b), and (c). 

(a) V{-oo) = Vi. 

Substituting (3.11) into (4.2), we obtain 

Sr{x, k)-^^a^ + ikaf + {ikfaf + ■■■ + {ikfa^ + 5"^, (4.3) 



a, 



R _ ^ 1™ „-W+V(x) (^a IN „R 1 







lim e-^+^(^) (f^ - 1) , a^=- lim e-^+^(^) (n > 1), (4.4) 



4 W-^-oo V u y "4 w-^-oo 



hm e-^+^(^)p^. (4.5) 
Prom (3.12) and the first equation of (4.4), we have 

o^{x) = -le-^ie^(^). (4.6a) 
2 

For n > 1, we obtain by substituting (3.13) into the second equation of (4.4). We 
can see from (3.15) that PiS.,<T„_i = if m < A. So m - A > 1 in (3.14). Taking the 
hmit W — >■ — oo of e~^f^ amounts to picking out the terms with m = A-|- 1. This gives 

E ^ P^.....-^''' ^2, ■ ■ ■ , a^-i]-.. e^(-). (4.66) 

{(Tl-'-CTn-l} 

The explicit expressions for the first few n are 



(x) = {3e-^^^ (-++]^-oo - (— +]-oo}^ 
alix) = {l2e-^^^ (-+++]^oo - 2e-^^ (-+-+]^oo - Ge"^^ (— ++]-oo}^ 
af{x) = {60e-^^^ (-++++]^-oo - Ge"^^^ (-++-+]^-oo - ISe"^^^ (-+-++]^-oo 

- 36e-2^i (— +++]^_oo + 2 (— +-+]^oo + 6 ( ++]^_oo}e^(^\ (4.7) 

where we have used the shorthand notation {—]^^, { — i-]-ooj etc for {—1]^^, 
(-1,+1]^^, etc. 

Obviously \a^\ < oo if V—Vi G Fj^_{. (This is the same as the condition for |f^| < oo 
discussed in the previous section.) So, equation (4.3) makes sense ifV — Vie F^}^. 

(b) y(-oo) = +00. 

In this case, we substitute (3.18) into (4.2). This leads us to consider the hmits 

6^ = 1 hm e-^+^W(f^-l), ^^: = ^ lim e'^+^W (n > 1), (4.8) 

Sf = \jm e-^+^(^)p^. (4.9) 

4 W^—oo 
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Prom (3.19) and (4.8), we have = 0. For n > 1, we can see that only the terms with 
A = in (3.20) survive in the hmit of (4.8). Since A is an odd number for even n, it 
follows that — for any even n. So, in this case we have the expression 

Srix, k)-^^ ikbf + iik)^f + (ikfbf + ■■■ + iikf''+'b%^, + 5^1^,, (4.10) 

which has only odd powers of k. Since = {ik)"'b^ + S'^, we have 6^ = for 
even n. The coefficients b^ are easily obtained as 

bni^) [-1' e^(^\ (4.11) 



2 

{o"l-'°"n-l> 
A=0 



where the sum is taken with the constraint A = Y17=i — 0- Using the shorthand 
notation [— ]^oo, etc for [— l]^oo> ^^^j explicitly write, for the first few n, 

bf{x) = I [-]^_ooe^^^\ bf{x) = - [— +]!l^e^(^), 

bf{x) = {6 [—++r-^ + 2 [— +-+]^_oo}e^^^\ 

b^{x) = -{72 [ +++r_^ + 36 [ +-++r-^ + 12 [ ++-+]^_oo 

+ 12 [— +— ++]"_oo + 4 [— +-+-+]^oo}e''^"^- (4.12) 
Note that (4.12) can be formally obtained from (4.7) by letting Vi — >■ +oo and 

( ' ] — OO ^ [ ] — oo • 

It is obvious that < cxo if e^^ G Fji_\, which is the same as the condition for 



|r^| < oo studied in the previous section. Therefore, (4.10) makes sense if e e F21J . 
(c) V{—oo) — —00. 

In this case, we cannot obtain the expansion of ,5^ - I by substituting (3.22) into (4.2), 
since limvi^_^_oo e~^(fo — 1) is not finite. Instead, the expansion of {Sr — |)"^ can be 
obtained in the same way as in case (b). As shown in appendix D, we have 

^tkbf + {tkfbf + {tk^bf + ■■■ + {tkf^n%^, + 6%^, , (4.13) 



where b^ and are the quantities obtained form b^ and by changing the sign of 
the potential V. Namely, 

bnix) = l E ^iV.....-i[+l'-^i'-^2,...,-a„_ir_ooe-^^^\ (4.14) 

{<Tl---a„_i} 
A=0 

or, explicitly, bf{x) = \ [+f^^e-^^''\ b^{x) = - [++-]^^e-^(^\ etc. From (4.13), the 
expansion oi S^. — \ is obtained as 

Sr{x, k)-\^ {ikVl\ + iki^ + {ikfi^ + • • • + iikf''^'l2M+, + C+i, (4.15) 

m?-bfb^], etc. (4.16) 



1 b^ 



Abf A{b^y ^ A{bff 
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Obviously 7^ = and = S^^-^ for even n. 

We know that \b^\ < 00 if G F,j:^_\. (This is obtained from the condition for 
< 00 by changing the sign of V.) Since 7^ is expressed in terms of with 



m 



< n + 2, wc can see that |7^| < 00 if e -^i+i) hence that (4.15) makes sense 



V ^ p(-) 



The expressions for Si can be obtained in the parallel way. We consider the three 
cases, where V{+oo) is finite (= V2), +00, or —00. The results are as follows: 

(a) If V{+oo) = V2, 

1 



Si{x, k) - - = a^ + ika\; + {ikya^ + ■■■ + {iky a^^ + 



a\{x) 



2 



= E ^-^^^ i^ifX^?...,.„_,e-^^^ K-u . . . , -2, a„ -1; 

{(71---(7„_l} 

(b) If V{+oo) = +00, 

Si{x, k)-^^ ikh\ + {ikfh\ + {ikfh\ + • • • + (i^)'^+^6^M+i + <^2Vi 

^ni^) -\ P^IU..,^.-. ...,<72, <7i, -l]r e^(^) 

{<'l--<^n-l} 
A=0 



00 g^(^)_ 



(4.17) 

(4.18a) 
(4.185) 

(4.19) 
(4.20) 





-1 


Si{x,k) - ^ 


= 4 



(c) If V{+oo) = -00, 

ikb\ + {ikfh\ + {ik)% + • • • + (^A;)2^+^6^M+l + <^2Vi] , (4.21) 

Si{x, k)-\^ {ik)-'i\ + ik^\ + (iA;)^3^ + • • • + {ikf'^'^^M+i + (4-22) 

~^n(^) = ^ E ^iiU,...,.„-J-^n-l,...,-<72,-<7i,+l]re-^(^). (4.23) 

{0-l---(T„_l} 

A=0 



5. Low-energy expansion of the Green function 



The expansion of the function 5" (equation (2.6)) is obtained by adding the expressions 
for Sr (equation (4.3), (4.10), or (4.15)) and Si ((4.17), (4.19), or (4.22)). We can derive 
the expansion of by substituting this into (2.7). We study the six cases listed in the 
introduction. 

Case (i): V{-oo) = V^, F(+oo) = V2. 

Adding (4.3) and (4.17) together, we obtain the expansion of 5" — 1 as 

S -1 = Sq + iksi + {ikfs2 + ■■■ + {ik)^SN + 5jv, (5.1) 
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s„ = a^ + a^, S^^Sf + S^^. (5.2) 
Substituting (5.1) into (2.7) we can derive the expansion 

^8(2;, y; k) = {ik)-^g-i + go + ikgi + {ikfg2 + ■■■ + {ik)^ g^ + Ajv, (5.3) 

1 ( / \ 5i(x) s\{y) \ . ^. 

f-i = ^ / , , , T' ^ ^1^' ^ ~ o ( \ ~ o ( \ 9-u etc., 5.4 

2^yso{x)so{y) V 2so(x) 2so(|/)y 

where we have defined 



X 



(ln{x,y) = - Sn-i{z)dz. (5.5) 

(The remainder term A^v will be discussed in the next section.) Prom (4.6), (4.18), 
(5.2), and (5.4), we obtain the explicit expressions for g_i and gi, 

^-lV(x)+Viy)]/2 

[V{x)+V{y)]/2 r 1 / N ^ 

90 = ^ | (e-v-.+e-v-.)2 ( (H-oo + + (H^-oo + H^) + | > (5-66) 

which are the same as the expressions derived in [17] by a different method. 
Case (ii): V{-oo) = Vi, V"(+oo) = +oo. 

In this case, the expansion of 5* — 1 has the same form as (5.1), with 

sn^a^ + bl = C + (5-7) 

Since — and — for even n, equations (5.7) read s„ = and 5„ = + 
for even n. The expansion of Gs, too, has the same form as (5.3) with (5.4). Substituting 
(4.6) and (4.20), we obtain 

^_^^g-[m+y(.)]/2gy,^ (5.8a) 

-[V{x)+V{y)]/2 

90 = — 2 — {e^^^ ( {-r-^ + + (-i^-oo + [-ir) + w.^- (s-s^) 

Formally, (5.8) can also be obtained from (5.6) by letting ¥2^00 and [— )^ — )■ [—]^- 
Case (iii): y(-oo) = Vi, F(+oo) = -00. 

When V{+oo) = —00, the expansion of 5" — 1, which is obtained from (4.3) and (4.22), 
contains a term of order 1/k: 

S -1^ (iky^s^i + so + iksi H h (ik)^SN + 5n, (5.9) 

5-1=7-1, 5n = a^ + 7n (^>0), 5„ = 5f + e- (5.10) 

When n is even, s„ = since 7^ = 0. Corresponding to (5.9), the expansion of Gs 
lacks the term of order 1/k which was present in (5.3): 

Gs(x, y; k)^go + ikgi + {ikfg2 + • • • + {ikf gw + Ajv, (5.11) 
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-exp[qo{x,y)] f Soix) So{y) \ , . 

9o = ^ / , , 9i = Qiix, y) - 7, TT - o TT 9o, etc., 5.12 

where Qq and qi are defined by (5.5). For n — —1 and 0, we have = 
and So{x) = -e^(^)-^V2. Note that = -e^W. Therefore, 



ex.p[qo{x,y)] = exp 



po 



(5.13) 



Also using qi{x,y) = — So{z)dz = |e ^4+]^; from (5.12) we obtain 

^ _e-^vi.)+viy)]/2 [^joo^ ^ _e-[v(x)+y(.)]/2g-vi [^j^ [^]oo_ (514) 

Case (iv): V{—oo) — +oo, \^(+cx3) = -\-oc. 

When both V{±oo) are +oo, the expansion of 5" — 1 has only odd powers of k: 

iksi + {ikfss + {ik)^S5 + ■■■ + {ikf^^^S2M+i + Wi, (5-15) 

sn = hl + hl, <5„ = 5f + e. (5.16) 

(Note that Sn = and (5„ = 5„_i for any even n.) The corresponding expression for Gs 
begins with the term of order l//c^, and has only even powers of k: 

Gs{x, y; k) = (ik)-^g-2 + Qo + {ikfg2 + {ikfgA + • • • + {ikf^g2M + A2M, (5.17) 

^-2 = ^ / , X , T' 30 = (l2{x, y) - ^—r^ - g-2, etc. 5.18) 



(In deriving (5.18), we choose the branch of the square root in (2.7) so that \/—k'^ = 
—ik.) For n = \ and n = 3, the first equation of (5.16) reads 

s.{x) = ie^(^) [-]-^, s,{x) = -e^(^) ( [— + [+— ) ■ (5.19) 

(Note that [-]^oo + [-]^ = [-]-oo-) Hence q2{x,y) = [-\-oo- Substituting these 

expressions into (5.18) yields 

9-2 = -e-[^(^)+^(-)l/^pi^, (5.20a) 
^ e-[v(a;)+y(j/)]/2 I [ +]-oo + [+ ]g° + ]^ + Wl I . (5.206) 

( [-]-oo) 2 

Case (v): V{—oo) = +00, y(+oo) = —00. 

In this case, too, the expansion of 5* — 1 has only odd powers of k, but now the series 
begins with the term of order 1/k: 

S-l = {ik)-h_^ + iksi + {ikfs3 + . . . + {iky^+h2M+i + Wi, (5.21) 
S_i=7-1, «n = 6^ + 7n (^>0), Sn^S'^ + Sl"^. (5.22) 
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(For even n, we have s„ = and 5„ = (^n-i-) Correspondingly, the expansion of Gs 
begins with the term of order k^: 

Gs{x, y; k) ^ go + {ikfg2 + {ik^gA + • • • + {ikf^g2M + Asm, (5.23) 

-exp[go(a:;,Z/)] / , . si{x) si{y) \ 

9o = ^ / , . =r^^ 92 = q2{x, y) - ^ - ^ c/o, etc. 5.24 

2Vs-i(x)s_i(y) V 2s_i(a:) 2s^i{y)J 

We have s_i(x) = e^(^)/(2[+]^) and Si(x) = e^(^){[-++]^/([+]~)2 + Since 
s_i is the same as in case (iii), obviously go is the same as (5.14), i.e., 

^0 = -e-[^(-)+^(^)]/2 [+]-. (5.25a) 

To calculate ^2, we use /; e^(^){[-++]r/([+]r)'}(i-2 = and 
integrate by parts. Substituting this and Si into the second equation of (5.24) gives 

92 = e-[^(^)+^(^)l/2{ ( i-t^ + [-]^_oo Wy + [-+]y) + 2 [-++1^ }■ (5-256) 
Case (vi): \^(— oo) = —oo, V{+oo) = — oo. 

In this case, the expansion of 5* — 1 has the same form as (5.21), where 

«n = 7^ + 7n, Sn^SZ^' + SZ^ (5.26) 
(As before, s„ = and 5„ = Sn-i for even n.) The expansion of Gs is given by the same 
expression as (5.23) with (5.24). Substituting s_i(x) = e^(^)[+]!f^/ (2 H^ooW?") and 

exp[qo{x,y)] = [+]loc[+]T / {[+]-oo[+]^) , we obtain 

_ _p-[V{x)+V{y)]/2 [+]-oo[+]^ 271 

^0- e ^^j^^ . (5.27) 

(The expression for g2 does not become simpler than (5.24) in this case.) 

Now we have derived the expansion of Gs for each of the six cases. To summarize, 
the expansion has the form of (5.3) (in cases (i) and (ii)), (5.11) (in case (iii)), (5.17) (in 
case (iv)), or (5.23) (in cases (v) and (vi)). These expressions make sense if and only if 
\gn\ < oo for any n < N {N = 2M for (5.17) and (5.23)). The remainder term is 
finite if all the gn are finite. As mentioned in the previous section, we know that 

\a^\ < oo if V-Vie Fiz\, |a^| < oo if V-V2e F^t\, 



b^\<oc if e-^eFi:!, |6^|<oo if e-^'eFiti 

.V ^ jpi-) I L| ^ -r V ^ M 



|7^|<oo if e^eF^;l, |7nl<oo if e ^S- (5-28) 

Prom (5.28), we can derive the sufficient conditions for \g]^\ < oo. These conditions are 
also sufficient for Ig^^] < oo (n < A?"). As a result, we find that 

in case (i), equation (5.3) makes sense iiV — ViE F^^ and V — V2 E F^\ 

in case (ii), equation (5.3) makes sense iiV — Vi E F^^ and e -^jv^\ 

in case (iii), equation (5.11) makes sense if F — Vi e F^}2 and e -^jv"\ 

in case (iv) equation (5.17) makes sense if e P2M+2 ^'^'^ ^ ^2M+2^ 

in case (v), equation (5.23) makes sense if e F2M-2 and e -^2^, 

in case (vi), equation (5.23) makes sense if e F21J and e F2^ . (5.29) 
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(The conditions involving F^^^ with n < are interpreted as automatically satisfied.) 

6. Behavior oi as k 

The expansion of Gs to order is meaningful as a low-energy expansion only if the 
remainder term satisfies — o{k^) as /c — >■ 0. In this section, we study the conditions 
for this to hold. Substituting the expansion of — 1 into (2.7), we can easily see that: 

in cases (i) and (ii), A„_i = o(A;"-^) if 5„ = o(A;") (n > 0), (6.1a) 

in cases (in), (v), and (vi), A„+i = o(A:"+^) if (5„ = o(A;") (n > -1), (6.16) 
in case (iv), A„_3 = o{k''-^) if 5n = o{¥') (n > 1). (6.1c) 

(We have used the fact that the integral in (2.7) and the limit A; — > are interchangeable, 
as can be easily shown.) By definition, 5„ is equal to (i) + 5,"^, (ii) + 5^^^, (iii) 
df^ + (iv) + ^n^, (v) + or (vi) df' + according to the six cases. So, 
the behavior of Ajv as A; — )■ can be known from the behavior of 5"^, 5^^, etc. 

The small-A; behavior of 5^^, 5^^, etc can be studied using (3.16) and (3.21). The 
detailed analysis is given in appendix E. The result is: 

5f = o{¥') if V-V^eF^-X 6^^ = oik'') if V-V2eFi% (6.2a) 
(5f = o(fc-) if e-^eFi:l, C = o(n if e-^eFill, (6.26) 

Sl^i = oik-~') if e^GF(-), Sti = oik^-') if e^GFW. (6.2c) 

(Here n > 0. For n = 0, the conditions involving Fj:^'^\ should be interpreted as 
automatically satisfied.) From (6.1) and (6.2), wc can conclude that: 

in case (i), A^v 
in case (n), Ajv 
in case (iii), Ajv 
in case (iv), A^v 
in case (v), A^v 
in case (vi), A^r 

The conditions in (6.3) are exactly the same as as the conditions in (5.29) (where 
N — 2M for cases (vi), (v), (vi)). Therefore, the expansion to order makes sense 
and is valid as an asymptotic expansion if these conditions are satisfied. 

The marginal cases for the conditions of (6.3) are V{z) ~ A + /3/\z\°' (where a, /3 
are constants, and A = Vi or V2) and V{z) ~ a log \z\ as z ^ —00 or +00. These cases 
correspond to Vs{z) ~ C/|zp+" and Vs{z) ~ /(/ + l)/|-2p with I = a/2, respectively. 
Let q; be a non- integer such that 0<n<Q;<n-|-l. Then, as shown in appendix E, 

5^^~CA;" (A;^0) if V(z) r^Vi + j^ (z ^ -00), (6.4a) 
5^^ ~ CA;" (k 0) if V(z) ~ a log \z\ (z -00), (6.4&) 



= o{k^) if V-VieFj^> and V-V2eF}^> {N > -1), 

= o{k^) if V-Vie Fj~^ and e"^ e Fj^^ {N > -1), 

^o{k^) if V-VieF^j^}^ and G {N>0), 

= o{k^) if e-^ e FjSri and e"^ e Fj^^^ {N > -2), 

^o{k^) if e-^eFj^}^ and e F^^^ {N>0), 

-o(fc^) if e^eF^-^ and e f}^^ (N>0). (6.3) 
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and similarly for 5^ and S'^. (Here C is a certain constant.) The behavior of Ajv for 
the marginal cases can be easily known from (6.4) (and the corresponding expressions 
for and 6^). For example, if — Vi G -Pat+i and V{z) ~ a log 2; as z ^ +00 with 
A^ + l<Q!<7V + 2, then Ajv ~ Ck"~^ as k Q (see example 5 of section 8). 

7. Schrodinger equation with a potential vanishing at a; — >■ ±00 

Suppose that Vs{x) is given, and that f{x) and V{x) are yet unknown. Let ipoi^) be a 
solution of (1.1) with A; = 0. Then a function / satisfying (1.4) is obtained from i/jq as 

f{x) = ^\ogM^)- (7.1) 

The Schrodinger equation (1.1) is equivalent to the Fokker-Planck equation (1.2) with 
(7.1), where is related to ip by (f){x) = iIjq{x)iIj{x). The Fokker-Planck potential is 
expressed in terms of ipo as V{x) = — 21og'?/'o(a;)- (In order to make V{x) finite for any 
finite x, the function ipo needs to satisfy ipo{x) > for any finite x.) A given Schrodinger 
equation can be thus transformed into a Fokker-Planck equation. 

In our formalism, Vs{x) need not be zero (or even finite) as x — > ±00. But here we 
study a particular feature of the case where Vs{x) tends to zero at both x — > -|-oo and 
X — >■ —00. For simplicity, we assume that there are no bound states. 

Now we assume Vs(±oo) = 0. Let iPq{x) and iPq{x) denote the solutions of (1.1) 
with k — 0, such that iPq{x) ^ 1 as x ^ -|-oo and ipo{x) 1 as x ^ —00. We define 

V±(x) = -2\og^^(x), /±(x) = ^log^^ix). (7.2) 

Obviously, V+{x) — )■ as a; — )■ +00 and V^{x) — )■ as x — )■ —00. If ipQ and ip^ are 
linearly independent, then fj^^f^. By using or /_ in (1.2) in place of /, we have 
two different Fokker-Planck equations equivalent to (1.1). (As a matter of fact, we can 
take any linear combination of ■^J" and ■0^, so there are an infinite number of equivalent 
Fokker-Planck equations.) If and %1)q are linearly dependent, then /+ = /_. In 
the conventional terminology of scattering theory, the cases /_|- 7^ /_ and /+ = /_ are 
referred to as "generic" and "exceptional", respectively. In the exceptional case, V-{x) 
(= V^{x) + constant) is finite at both x — )■ +00 and x — )■ —00. Thus, the exceptional 
case corresponds to case (i) in our classification in section 5. In the generic case, on the 
other hand, V-{x) and V+(a;) tend to —00 as x ^ +00 and x —00, respectively. So, 
the generic case is included in our case (iii). 

There is no particular difficulty in dealing with the exceptional case by our method; 
we can directly use the results of section 5 for case (i). On the contrary, special care 
is needed for the generic case. In the generic case, iPq{x) grows linearly, and so V-{x) 
diverges logarithmically, as x +00. We can sec that e^- ^ since e behaves 
like 1/x'^ as X — >■ -|-oo. The criterion (6.3) for case (iii) indicates that the expansion to 
order may not be valid for > 1 if we use V- in place of V. (The situation is the 
same for V-^, since ^ Fi~'* ■) Fortunately, we can avoid this difficulty by using both 
V+ and V-, as explained below. The idea is to use V- for Sr, and V+ for Si. 
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The relations (2.10) hold for both V+ and V^, so 

Rr _ c+ I f+ _ c- I /- _ q+ f+ _ Q- f- /7 n\ 



1 + Rf ' 2ik ' 2ik' 1 + Rf ^ 2ik ' 2ik' 
Here 5"^ and denote 5",. and 5"; with V± in place of V. Prom (2.11) and (7.3) we have 

S{x, k) = ^ [f.{x) - /+(x)] + S;{x, k) + St{x, k). (7.4) 

The expansion of 5* — 1 takes the form of (5.9), where, instead of (5.10), 

s_i = l[/_(x)-/+(x)], s, = a^ + a^+ (n>0), 5^ = 5f- + 5f+ . (7.5) 

Here and are defined by (4.6) and (4.18) with V- and respectively, in place 
of V (and similarly for 61^~ and 5f^). In particular, sq = -|(e^-(^) + e^+(^)). (Note 
that Vi = y_(-oo) = and V2 = V+{+oo) = 0.) Substituting (7.5) into (5.12) gives 

_ -exp{[V_ix)-V^ix)-V4y) + V4y)]/A} 

Vif-i^) - U{x)] [f-{y)-Uiy)] 



1 

f 1 = o 



V- (x) j^^V+(x) g {y) _|_ gy+ (y) 



A f-{x)-f+[x) f-{y)-f+{y). 



go. (7.6&) 



The higher-order coefficients can be calculated without any difficulty by using (7.5). 
Now (5.11) makes sense, and = o{k^) as A; ^ 0, if e F^}^ and V+ e F^^}^. 

8. Examples 

To demonstrate the calculation of the expansion, let us consider some simple potentials 
for which the exact form of the Green function is available. (For examples 1-4, the 
expressions for the exact Green function can be found in section 11 of [4]. Note that 
G = 2ikGs in [4].) In all the graphs, k is taken to be a real number {k >0). 

Example 1. V{z) = z^, Vs{z) = z"^ -I. 

The first example is a parabolic potential. The corresponding Vs is also parabolic. From 
(5.16), (4.11), and (4.20), we obtain the first two coefficients of (5.15) as 



, , 0r 2 TT 



/'OO n 

I {eric zY dz -\- I {ericzYdz 



(8.1) 



and form (5.5) we have q2{x,y) = (7r/4)(erfi|/ — erfix). (Here erfcz = 1 — erf 2; = 
(2/v^) j^e-"" dw, and erfiz = -iexl{iz).) The expansion of the Green function has 
the form of (5.17), where gr_2 and qq are obtained by substituting the above expressions 
into (5.18). The approximation to this order, Gs — {ik)~'^g-2 + go, is plotted in figure 1 
along with the exact value. (The exact expression of the Green function is given by 
equation (11.6) of [4]. Note that this Gs takes a real value when k is real.) Higher-order 
coefficients of the expansion of Gs can be obtained in the same way. Although we omit 
here the expressions for S5 and g2, the result of the calculation up to order k'^ is also 
shown in figure 1(a). 



Low-energy asymptotic expansion of the Green function 



18 




Figure 1. The Green function Gs{x,y;k) for tlie potential V{z) = (example 1), 
(a) plotted as a function of k, with x = 1.2, y — 1; (b) plotted as a function of x, with 
k = 1, y = (here Gs is real). In all the graphs, solid lines are the exact values. Here 
the dashed lines show Em=-i(«fc)^".92m- {N = 2M = -2, 0, 2 in (a), and = -2, 
in (b).) In (a), the curve labeled as 'W — 2 (modified)" is the plot of (8.2). 



The exact Gs{k) has poles at k = ±y/2n [n = 0, 1,2, . . .), corresponding to the 
eigenvalues of the Schrodinger operator (/c^ = 2n). From the information of go and g2, 
we can approximately reproduce the poles nearest to the origin as 

'^s - TTT^9-2 + , ^ ^ (8.2) 

This is a better approximation than Gs — (^^)^^5'-2 + ^70 + {ik)^g2 (see figure 1(a)). 

Example 2. V{z) = 2 log cosh z, Vs{z) = 1 — 2sech^2;. 

In this example, V{z) diverges to +oo and Vs{z) tends to 1 as ^ — > ±oo. As in the 
previous example, the expansion of Gs has the form of (5.17). From (5.16), (4.11) and 
(4.20), we can easily calculate 

Si(x) = cosh^x, S3(x) = — - cosh 2x cosh^ X, (8.3) 

and q2{x,y) = |(?/ — x) + |(sinh2|/ — sinh2x). By substituting them into (5.18), we 
obtain g-2 and go. The higher-order coefficients can be similarly calculated. The results 
are shown in figure 2. (For the expression of the exact Gs, see equation (11.15) of [4].) 

This Gs has branch point singularities at A; = ±1, and ImGs = for |A;| < 1. The 
series ^m=-i(^^)^™5'2m convergent for \k\ < 1. For \k\ > 1, we can use the high-energy 
expansion (discussed in [4]) to calculate Gs with very good precision (see figure 2). 

Example 3. V{z) = e^ Vs{z) = \e^^ - \e\ 

This exponential potential belongs to case (ii) of section 5. The expansions of S and Gs 
have the form of (5.1) and (5.3), respectively. From (5.7), (4.6), and (4.20) we have 

s^{z) = -iexp(e^), s,{z) = exp(e^)[Ei(-e^) + 2Shi(e^)], 

S2(z) = -2exp(e") / exp(e'")Shi(e"') du;, (8.4) 
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exact 



ReGs ImGs N=-2 




Figure 2. The real and imaginary parts of Gs{x,y;k) for the potential V{z) = 
2 log cosh 2: (example 2), plotted as functions of k, with x ^ 2 and y — 0- The dashed 
lines show E^^„i(ifc)2™g2m. {N = 2M = -2, 0, 2.) The dotted lines are the result 
of the high-energy approximation (see [4] ) obtained by expanding log Gg in powers of 
1/k to order (For fc > 1, the curves of the high-energy approximation almost 

coincide with the exact curves.) 



and so on, where Ei and Shi denote the exponential integral function and the hyperbolic 
integral function, respectively. (Ei(z) = — J^(l/t)e~* c/t, Shi(z) = J^{l/t)smh.tdt.) 
We obtain g^i, go, etc by substituting (8.4) into (5.4), etc. The results of the calculation 
to order k"^ are shown in figure 3. (See equation (11.12) of [4] for the exact form of Gs-) 
With fixed y and k, the function Gs{x,y; k) rapidly falls off to zero as a; — ?■ +00, 
and oscillates as x — )■ —00. (Recall that Gs{x,y) = Gs{y,x). Since we have been 
assuming y < x, the Green function for x < y is Gs{y, x) with our expressions.) As can 
be seen from figure 3(b), the truncated series {ik)~^g_i + go + ■ ■ ■ + {ik)^ gj^ gives a good 
approximation of G^ as a function of x for x > y. (In figure 3(b), the approximation with 
N = 2 almost coincides with the exact value for x > y.) However, this approximation 
is not effective when x < y and y — x is large. To cope with the oscillatory behavior of 
Gs as X — )■ —00, it is better to truncate the expansion of logGs (rather than Gs itself), 
and then exponentiate it. Namely, 

Gs - {ikY^g-i exp [ikpi + {ikfp2 + ■■■ + (^fc)^+ Wi] , (8.5) 

where {pi, . . . ,pn+i} can be expressed in terms of . . . , (^tv} as pi = go/g^i, 

P2 = {gi/ g-i) — \{go/ g-iY , etc. As shown in figure 3(c), this gives a good approximation 
in the region where — x is large. 




V.iz 



- zY' + \{i - zY"' (.<o) 



yiT^+i, 1 ^(i + ^)-i-i(i + z)-3/2 (^>o). 



This example belongs to case (v) of section 5. Here V{z) slowly diverges to =Foo as 
z — )■ ±00, and Vs(2;) tends to zero like \z\~^ . In this case, it is easier to use (5.25a) and 
(5.25b) directly for the calculation of go and g2. Assuming that y < < x, we obtain 



^0 



-2exp(l-^-^)(l + ^). (8.6a) 
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exact N=0 

N=2 



Figure 3. The real and imaginary parts of Gs{x,y;k) for the potential V{z) = 
(example 3), (a) plotted as functions of k, with x — 0.5 and y = 0; (b) plotted as 
functions of x, with k = 0.4 and y = 0. The dashed lines show where 
= 2 in (a) and = 0, 2 in (b). (Since k is real, = and N — 2 are the same 
as = — 1 and A^ = 1, respectively, for the imaginary part.) In (c), the same graphs 
as (b) {Gs{x) with k — 0.4, y = 0) are drawn with a larger scope. The dashed lines 
in (c) are the plots of (8.5) with A^ = and 2. (They are plotted only for x < 0.) 



5(2 = - exp 1 



118 + 37x + 2x^ -3?/ 



+ (94 + 11a; - 3y)y/l + x + 2(1 - y) (l + Vl + x ) y^l-y . {8.6b) 

As can be seen from figure 4, equation (5.23) with (8.6) gives the correct asymptotic 
expansion of the Green function. (The exact Green function is given by equation (11.25) 
of [4] with the replacement (x, y) — ?■ (— —x). The V{z) in this example is the same as 
V{~z) in example 7 of [4].) The series {ik)~'^g_2 + 90 + (^^)^5'2 + ■ ■ ■ takes a real value 
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Re Gs Im Gs 




k 



Figure 4. The real and imaginary parts of Gs{x, y; k) for the potential of example 4, 
plotted as functions of k, with x — 1 and y = —0.5. The dashed line is the plot of 
go + (ifc)^(72- The falloff of |ImGs| as fc — )■ is faster than any power of k. 



ReGs ImGs 




k 



Figure 5. The real and imaginary parts of Gs{x,y; k) for the potential of example 5 
with a — 3/2, plotted as functions of k, with a; = 1.5 and y = 0.8. The difference 
between the solid and dashed lines is Aq. Here Aq behaves like fc"^^ — y/k as fc — > 0. 



when k is real. Although the exact Gs is not real, the imaginary part of it approaches 
zero faster than any power of as A; — )■ (see figure 4). Since Gs{k) is essentially 
singular at /c = 0, the series X]^^L-i(^^)^"fl'2n is asymptotic but divergent. 

Examples. V{z) = ae{z - 1) log z, Vs{z) = ^^^^e{z ~ l){l/z^) - 16{z ~ 1). 

Here a is a positive constant, and 6 denotes the heaviside step function. {V{z) = for 
z < 1.) Since f{z) is discontinuous at z = 1, the Schrodinger potential contains a delta 
function at 2; = 1. Assume that y < 1 < x. From (5.8) we obtain 

9-1 = 90 = [l~y + ^dz^ . (8.7) 

Obviously, go is finite if a > 1 (e~^ G Fq~^^) and infinite if a < 1 (e~^ ^ -^o^^)- '^^^ 
Green function for this potential can be exactly obtained as 



where is the Bessel function. Since J^{z) behaves like z'^ as z ^ for non- integer u, 
we can see that the asymptotic expansion of (8.8) contains a term proportional to 
IfA^ + l<a<A^ + 2, then ~ C/c"~^, as explained in section 6 (see figure 5). 
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ReGs ImGs 




Figure 6. The real and imaginary parts of Gs{x,y; k) for the potential of example 6 
with a = 1, plotted as functions of fc, with x = 0.5 and y — —0.5. 



Example 6. 

'a^ (1^1 < 1) 



i\z\>l), 



where a is a constant. This example is to demonstrate the method discussed in section 7. 
The zero-energy wave functions ^jj^ for this Vs are 

r 1 (C1-C2X {x< -1) 

i)Q{x) = < cosh[a(a; + 1)] i^oi^) = { cosh[a(x - 1)] (-1 < a; < 1) (8.9) 

[Ci + C2x, [1 (Kx), 

where Ci = cosh(2a) — asinh(2a), C2 = asinh(2a). We obtain V± and f± from (7.2). 
Let us consider the case —l<y<x<l. Substituting the expressions for V± and f± 
into (7.6a) and (7.6b) yields 

cosh[a(x - 1)] cosh[a{y + 1)] /o 1 n ^ 

90 = ^TTT^ ' 8.10a 

a smh(2aj 



I tanh[a(x + 1)] + tanh[a(x — 1)] — tanh[a{y + 1)] — tanh[a(?/ — 1)] 



9o 

'' = 2-a\ 

1 J cosh[a(x — 1)] cosh[a(x + 1)] cosh.[a{y — 1)] cosh[a{y + 1)] 
sinh(2a) \ cosh[a(x + 1)] cosh[a(x — 1)] cosh.[a{y + 1)] cosh.[a{y — 1)] 

(8.106) 

On the other hand, the exact Green function is obtained by a standard method as 

^ ^ [{p- ik)eP^^-''^ + {p + iA;)e-P(^-^)] [{p + ik)e-P^^+y^ + {p - ik)eP^^+y^ 

^ ' ^' —4p [{p^ — k'^) sinh(2j9) — 2ipk cosh(2p)] ' 

p = Va^- k^. (8.11) 

It is not difficult to check that (8.10a) and (8.10b) are the correct coefficients of the 
expansion. The higher-order coefficients can be calculated by using (7.5) (see figure 6). 
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9. Summciry and remcirks 

The asymptotic expansion of Gs in powers of k is obtained by substituting the expansion 
of S into (2.7). Since S — Sr + Si, we can treat Sr and Si separately. When V{—oo) is 
finite or +oo, the coefficients of the expansion of Sr can be obtained in a simple form for 
arbitrary order of k (equations (4.3), (4.6), (4.10), and (4.11)). When V{—oo) — — oo, 
the expansion of {Sr — |)~^, instead of Sr, takes a simple form (equations (4.13) and 
(4.14)). In a parallel way, corresponding expressions are obtained for Si (equations 
(4.17)-(4.23)). The behavior of the remainder term, and the validity of the asymptotic 
expansion, can be studied by using the expressions (3.16) and (3.21). The result is given 
by (6.3). For the "generic" case with Vs(x) such that Vs(+oo) = Vs(— oo) = 0, we need 
to use a modified method explained in section 7. In this case, too, the expansion of S 
to arbitrary order can be obtained in a simple form (equations (7.5)). 

A different method for the low-energy expansion of Gs is discussed in [17]. The 
method of [17], which does not use the reflection coefficients, is more direct than the 
method of the present paper. However, the derivation of the expansion in [17] is formal, 
and does not provide a way to estimate the remainder term. Practically, the methods 
of [17] and the present paper complement each other. It is a future problem to unify 
these two methods. 



Appendix A. Proof of (2.10) 

When a — — oo, the definitions for / and Vs (equations (2.2) and (2.8)) read f{x) — 
f{x)9{b — x) and Vs{x) — Vs{x)9{b — x), where 9 is the Heaviside step function. So, 

Pix) + fix) = Vsix) - f{h)5{x - h). (A.l) 

Namely, the Schrodingcr potential corresponding to / (equation (1.4) with f ^ f) 
differs from Vs by —f{h)5{x — h). It is an elementary exercise in quantum mechanics 
to show that the transmission coefficient for the Schrodinger equation with the delta 
function potential —f{b)S{x — b) is To = 2ik/[2ik+f{b)], and that both the right and left 
reflection coefficients are equal to Rq = —f{b)/[2ik + f{b)]. The reflection coefficient 
for the Fokker-PIanck equation includes the multiple reflections caused by this delta 
function. We can take the sum of these multiple reflections as 

Rr{b, -oo) = Ro + T^R'rib, -oo) J2 [RoR'rib, -oo)]" = Ro + . p ps;, ^ 

^ 1 - RoR^{b, -oo) 

^ -f(b) + [2tk-f{b)]Rf{b,-oo) 

2ik + f{b) + f{b)RS,{b,-oo) • ^ ■ ' 

Prom (A.2) we obtain Rr/{1 + Rr) = Rl/{1 + R^) - f/{2ik). The second equation of 
(2.10) can be proved in the same way. 
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Appendix B. Calculation of and 

From (3.12) and (3.4) we have 

a ,W-Vx , W -Viz) 
rl^i{z,W) =-tanh ^-^ + tanh 

smh[yi - V(z)\ ( W -Vi V^- Viz) \ ^ , 

= — -j-^ Ht tanh ^ + tanh^ ^ . B.l 

sinh[W^ - V{z)] V 2 2 y ^ ' 

Substituting (B.l) into = C{fQ + ^) gives 

fl(x, W) = 2 J" dz-^ {siiihUy - V(z)]lf'„ + mz. »'))} 

= 2£d.sinh[V. - V(z)]^t^^'^^ = ^^^^ (B.2) 
We can express (cosh^ ^~^^ ) ^ as an infinite series in powers of and write 

oo 

t{x,W) = 2[^(-l)™+ime™(^-^i)]e^^ (-l]-oo- (B-3) 



oo 
m=l 

Let us define the operators 



^(2) _ / 1 , ^ \ ^(2) _ -W ( - ^ 



Then equation (3.5) can be written as 

L gix, W)^Y1 r ^"""^'^J-a 9i^, W) dz, (B.5) 

where j!:} stands for J'!i^ and for cr = +1 and a — —1, respectively. Substituting 
(B.3) into = r^-iff , and using /^^ dze"^^'^ {■ ■ -Jl^ = (•••, (j]-oo> we obtain (3.13) 
with 

oo 

D.,,.,,...,^„_,{W) = 2e^^ • • • J^lj^l 5^(-l)-+^me-(^-^^). (B.6) 

m=l 

Since jlje"^ = (-o-)(n - (7)e("-'^)^, it is easy to see that 

7-(2) . . . -7(2) pmW _ p(m) p(m-A)W /o yN 

with Pi\^l.,a„-i defined by (3.15). Hence we have (3.14). Expressions of and 
without the infinite sum over m can be obtained by carrying out the calculation using 
the last expression of (B.2) instead of (B.3). 

For the case V{-oo) = +oo, we have sinh[l^ - K(z)][f^ + ^{z, W)] = e^-'^W - 1. 
Substituting this into — £(fQ + ^) (see the first hne of (B.2)), we have 

r1 = 2 £^ dz^ (e^-^W - 1) = 2e^[-l]^_. (B.8) 

Hence we obtain 

f'„{x,W) = 2 Yl ill-,---^-^2.e^[-l,<7i,a2,...,a„_ir_. (B.9) 
Equation (3.20) follows from (B.9) and (B.7). 
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Appendix C. Finiteness of and 

When V{-oo) = Vi is finite, we have e±^(^) < C and | smh[Vi - V{z)]\ < C\V{z) - Vi\ 
for — oo < z < X. (Here and hereafter C denotes a constant which may not necessarily 
be the same everywhere.) Using these inequahties in (3.2a), we find 

\{-l,ai,a2,...,an-it^\ <C f |V^(^) - ^il 1^1""^ < oo (C.l) 

J — oo 

if V -Vie From (C.l) it follows that |f^| < oo if \/ - G F^j. 

Next, let us consider equation (3.20) for the case V{—oo) = +oo. Let us assume 
that Pau...,an-i 7^ 0, and let {pi,P2, ■ ■ ■ ,Pm} be a subset of {1, 2, . . . , n — 1} such that 
(Tp. = +1 for each i. Then there exist {^i, q2, ■ ■ ■ , Qm} £ {1; 2, . . . , n — 1} such that 

(Tq. = —1 and Qi < Pi for each i. (Otherwise Pal]...,crn-i = 0, as can be easily seen from 
(3.15).) If -Zp^ is sufficiently large, V{zq.) > V{zp.). So, exp[-V (zg^) + V (zp.)] < C for 
any Zg. < Zp. < x. We also have exp[— V^(z)] < C for any z < x. Therefore, 

|^i:L.„-i[-l'^i'^2,...,an-if_oo| r e-^(^)|zr^dz<oo (C.2) 

^— oo 

if e"^ e j. Prom (C.2) it is obvious that |f^| < oo if e"^ e F^^zl- 



Appendix D. Proof of (4.13) 

We use the rotation of coordinate axes discussed in section VI of [16]. From equations 
(D.l) and (D.2) of [16], we have Rr,-n-iT/2 = — (-Rr,o,-7r/2)"^, where Rrfi,e' is defined by 
(10.2b) of [16]. When V{—oo) — — oo, the Rr{x, —oo) corresponds to it!r,o,-7r/2- (Note 
that 9' — — 7r/2 corresponds to V{—oo) — — oo.) Therefore, 

1 \ ~ 1 ~ Rr,0,-n/2 ( c, l\ — -Rr,7r,-7r/2 



-4Lg,-- =2 -r,.,-.,. ^ [Sr--] =2 "------/^ (D.l) 

Thus, the expression for [Sr — has the same form as the expression for —4(5'^ — |) 
with 9 — instead of ^ = 0. As can be seen from (6.1) and (6.3a) of [16], the rotation 
with angle 9 — amounts to changing the sign of V and k. Hence we obtain (4.13). 



Appendix E. Estimation of the remainder term 

Here we consider (5^^ and (5^^, and prove the first halves of (6.2a) and (6.2b). The 
conditions for 5'^ and 5^^ can be derived in the same way. Equations (6.2c) easily 
follow from (6.2b) by noting that dfl^ = o(A;"-i) if 5%^ = o(A;"+i), and that is 
obtained from 5^ by the replacement V — )■ —V. 

Let us first note that 5'^ = o{k"') and 5^ — o{k^) hold, respectively, if — o{k'^) 
and = o{k"') hold for W — V{x). We can see this by substituting the expansion 
of Rr{x, —oo;k) with the remainder term (which is obtained by setting W = V{x) in 
(3.11) or (3.18)) into the first equation of (2.5). Here we show, more generally, that 
= o(/c") and = o{k"') hold for any W under the conditions stated in (6.2). 
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For k = 0, the scattering coefficients can be exactly obtained. (Altliougli cqiiations 
(2.3) liold only for k ^ 0, we can define the scattering coefficients for /c = by taking 
the hmit k 0.) We have [16] 

W - V(z) 

f{x, z;W; k ^ 0) ^sech (E.la) 

W -Viz) 

Ri{x, z;W;k^O)^ -Rr{x, z;W;k^O)^ tanh {E.lb) 

In this appendix, we make use of (E.l) together with the asymptotic expressions of f 
and Ri given in appendix F. 

Now let us derive (6.2a). For the case V{—oo) — we write (3.16) as 

ikf^^ E r dzA{z,h){-\,a^,...,a^_^Y_^e"-''^'^\ (E.2) 

|(Tl,...,<Tjv} 

where A{z, k) denotes the part containing f , Ri, and the sum over m. (We omit to write 
the dependence onW.) Equations (E.l) yield 

f^{x,z;k = 0) _^ 1 + Ri{x,z;k = 0) _ ^y^,_y^^^ 



P% 



1-Rf{x,z]k^0) ' 1- Ri{x,z;k ^0) ' ^ 

From equations (F.23) and (F.24) of appendix F, and from the fact that C3, C4, 
and 9 in these equations remain finite as /c — > 0, we find that |f^/(l — ^f)| < C, 
|(1 + Ri)/{1 — Ri) \ < C for —00 < z <x and < k^ with some k^. (Here, too, we let 
C denote a constant which may not necessarily be the same at each appearance.) Using 
these inequalities in (3.16), we can see that 1^4(2;, A;) | < C. (The infinite sum in (3.16) 
does not cause any problems. See the comment above (3.17).) In the same way as in 
appendix C, it can be shown that 



J — c 



dz|(-l,ai,...,ajv-i]looe'^-^(^)| <oo (E.4) 



\iV-Vi G F^j \ Since \A{z, k)\ < C, inequality (E.4) means that the absolute value of 
the integrand on the right-hand side of (E.2) is dominated by a /c-independent function 
of z which is intcgrable in the interval (— oo,a;). Therefore, if — Vi G -Fj^ we can 
interchange the order of the limit k ^ and the integral in (E.2) to obtain 



&o7a^= ^ / dzA{z,k^O){-l,a^,...,a^^^r_^e'^-''^^\ (E.5) 

|o-l,...,(Tiv) 

Substituting (E.3) into (3.16), and comparing it with (3.13) and (3.14), we find that the 
right-hand side of (E.5) is equal to f%_^^. Replacing ^ n — 1, we can write (E.5) as 



Since Pn_i = (^fc)"^" + /o", from (E.6) it follows that = o{k'^), and hence = o{k"'), 
as /c — > if V — Vi G Fji_\. (The proof for n = can be done by replacing with 
Rr, and using the integral representation of Rr{x, —00) (equation (3.16) of [3]).) 
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If V{z) — Vi ^ (^/\A'^ (A^ < a < iV + 1) as 2; — >■ —00, we can say more about the 



C\z 



N-a 



as z 



—00. 



behavior of 5'^ as /c — )> 0. In this case (—1, (7i, . . . , (T7v-i]1o< 
From (F.21) and (F.22) of appendix F, we can see that the leading contribution to the 
integral of (E.2) has the form 



r dzA{z, k) (-1, (71, ... , (7jv-i]iooe'^-^(^) = c r 

J —00 J —(. 



dz h{e 



ikz\ 



\N-a 



+ 



(E.7) 



where h is some function. The integral on the right-hand side is convergent if /c 7^ 0. 
This integral behaves like as /c ^ 0, as can be seen by changing the integral 

variable from z to Z = kz. So we have pN ~ Ck°', and hence ~ Ck°', as A; — )■ 0. 
Let us proceed to (6.2b). For the case V{—oo) — +00, we write (3.21) as 



(E.8) 



r 1 J —00 

{(71,...,(7jv} 

where B{z, k) is the part containing f and Ri. Let us temporarily assume that the order 
of the limit /c — )■ and the integral in (E.8) can be interchanged. Then, 



Pn 



lim , , , -, 
fc^o {ik)^+^ 



E 

{<Ti,...,(TAr} 



f dzB{z, 
J —00 



k^O) [-1,(71,..., (7jv-i]looe 



(E.9) 



Substituting (E.3) into (3.21) and comparing it with (3.20), we can see that the right- 
hand side of (E.9) is equal to r%^i. If G then f^+i is finite (see appendix C), 
and so (E.9) makes sense. With the replacement N ^ n — 1, equation (E.9) reads 



Pn-l 



lim 

fc^o [ik) 



if 



(-) 



n-l- 



(E.IO) 



Since p^_]^ = {ik)'^rl^ + p^^, it follows from (E.IO) that — o{k^), and hence S'^ — o{k'^), 



as A; if e"^ G F 



n— 1 • 



Now wc have only to justify the interchanging of the limit /c — )■ and the integral 
in (E.8). This is easy if f{—oo) ^ 0. When f{—oo) 7^ 0, the behavior of f{x,z) as 
z — > —00 is given by either (F.9) or (F.16). These equations hold for A; = 0, too, since 

f(x,z;fc = 0) = 2e^/V^(^)/2^1 + o(l)] (^ ^ -00), (E.ll) 

as can be seen from (E.la). The expressions (F.9) and (F.16) continuously approach 
(E.ll) as /c — )• 0. Since the quantities rj{z,k) and 9{z,k) in (F.9) and (F.16) are o{\z\) 
as ^ ^ —00, we have, for |A;| < A;o with some ko, 

\f\x,z;k) \ <Ce^p[-V{z) + C'z], (E.12) 

where C" is a constant which can be chosen arbitrarily small. Considering the behavior 
of Ri given either by (F.IO) or (F.17), we see that 



1 



1-Rf \1-R, 



l + Ri 



l-A 



1 



\l-Ri 



l + Ri 



l-Ri 



<C, 



since A < 0. From (E.13) and (3.21), we find 
\B{z,k)\ < Cf\x,z;k) e(^-^)^W. 



(E.13) 



(E.14) 
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In the same way as in appendix C, it can be shown that* 

|[-l,ai,...,a^_i]i^e'^-^(^)e(i-^)^(^)| <Ce^(^) T e-^'^^^lwf-Uw < Clz]"^ . (E.15) 

•/ — oo 

These inequahties hold as long as e F^}^. Prom (E.12), (E.14), and (E.15) we have 

\B{z,k)[-l,ai,...,aN-iY^ooe''''^^"^\ < C\zf cxp[-Viz) + C'z]. (E.16) 

Since /(— oo) 7^ 0, in this case V{z) tends to +cxo linearly or faster as 2; — )■ —00. 
Therefore, the right-hand side of (E.16) is integrable in the interval {—oo,x). The 
absolute value of the integrand of (E.8) is thus dominated by a A;- independent integrable 
function, and this justifies the interchanging of the limit and the integral. 

When V{—oo) — +00 and /(— 00) = (i.e., when V{z) grows slower than hnearly), 
we cannot find a /^-independent integrable function of z that dominates |r^(x, z; k)\ as 
in (E.12). In this case, the behavior of f{x, z) as z —00 is given by (F.23). However 
small k may be, f{x, z; k) is considerably different from f(a;, z\ 0) when —z is large, since 
(F.23) is not compatible with (E.ll). The crossover of the two different behaviors takes 
place in the region where \k\ ~ |/(^)|- (This can be known by studying the small-A; 
expansion off.) Let us define z^ by \k\ = \f{zk)\- (Such Zk is uniquely determined when 
k is sufficiently small, since we are assuming that f{z) is asymptotically monotone.) 
Roughly speaking, f{x, z; k) ~ f{x, z; 0) for z > Zk when k is sufficiently small. For 
z < Zk, we need to use (F.23). The factor C3 in (F.23) is of the order of e"^^^*)/^, since 
f(x, Zk, k) is of the same order as f{x, z^, 0) ~ Ce~^*^^'=)/^. We can write 

f(x, z; k) = Ce-^(^'=)/^ exp [-ik{z - Zk) + ie(zk, z, k)] [1 + o(l)], (E.17) 

where 9 is defined by (F.20). We divide the integral in (E.8) as /^^ = /f^ + The 
part can be treated in the same way as in the case /(— 00) 7^ 0. An inequality 
analogous to (E.12) holds for z > Zk, and it can be shown that limfc_^o /^^ = limjk->.o- 
Let us study the part /^^. Equation (E.17) givesji 

[x, z; k)\z\'' dz:^ Ce-^^"*) /A;^+^ (E.18) 



f2 



if we neglect the 9{zk, z, k) and the o(l) part of (E.17). (By using (F.18)-(F.21), it can 
be shown that the contributions from 9 and the o(l) part in (E.17) are indeed negligible 
in the limit k 0.) Using (E.15), (E.17), and (F.24), we can estimate the part of 
(E.8). This is essentially the same as (E.18), and we have 



Biz, k) [-1, a,,..., a^_i]l^e'^-^(^) dz 

-00 



< Ce-^^"-^ /k^'+K (E.19) 



(Also see (E.7) and the explanation below it.) If V{z) ~ (A^ + 1) log \z\ as z —00, 
then l^fcl ~ C/\k\ as k ^ 0. Since e'^^'^^ ~ C\zk\-^^+^^ ~ C\k\^+\ the right-hand side 
of (E.19) approaches a finite value as A; ^ 0. If e then V{z) grows faster 

* In appendix F of [3], it is assumed that the quantity on the left-hand side of (E.15) tends to a finite 
value as 2; — >■ —00, but this is wrong. 

1 When Imfc = 0, it is necessary to replace k hy k + ie and let e 4- after evaluating the integral. 
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than (A^ + 1) log \z\ a.s z ^ — oo, and so the right-hand side of (E.19) vanishes in the 
hmit k ^0. Therefore, the part /^^ of (E.8) is neghgible as A; ^ if e"^ G Fj^\ and, 
since hmfc_^o XT ~ /-oo ^^^k^o, this means that hmfc_^o and can be interchanged. 

If V{z) ~ alogl^l (A^ < q; < + 1) as 2; — 7> oo, the inequahties in (E.15) can 
be replaced by and (E.19) gives E ^(-2' ^) " " " ~ Ck"'^'^ {k 0). We 
can also see that B{z,k) ■■■dz^ C jf^\z\~''+^dz ~ C|^fe|-"+^+i ~ Ck"^-^'^ as 

A; — > 0. (This is obtained by substituting (E.3) into the integrand.) Therefore, from 
(E.8) we obtain p% ~ Ck'^, and hence 5^ ~ CA;", as A; 0. 



Appendix F. Asymptotic behavior of t{x,z) and Ri{x,z) as z —> —00 

In this appendix, we study the asymptotic forms of f{x, z; W; k) and Ri{x, z; W; k) as 
z —00 for Im/c > 0, /c 7^ 0. Details of the derivation are omitted, but let us only 
mention that equations (F.5), (F.14), (F.15), (F.21), and (F.22), which are the basic 
expressions, are all derived by using equations (3.5)-(3.8) of [16]. We need to consider 
the three cases, (1) /(— 00) = ±00, (2) /(— 00) = c 7^ 0, and (3) /(— 00) = 0. 

(1) /(-oo) = ±oo. 

Let us consider the Schrodinger equations 

-^^^(^) + [Tf'iz) + fiz)] i^^z) = fcV(^). (F.l) 

The equation for ip' is identical with (1.1), and the equation for is the Schrodinger 
equation corresponding to the inverted Fokker-Planck potential —V. We set 



i/^^iz) = exp 



1 

^Viz) + / p'^{w,k) dw 



(F.2) 



and substitute into (F.l). This gives the nonlinear differential equations 

-^p^iz, k) T 2f{z)p^{z, k) + [p±(z, k)f = -k^ (F.3) 
It is easy to see that these equations have solutions satisfying the asymptotic conditions 

p±(^,A;) = ±2^[l + o(l)] as z ^ -00. (F.4) 

We can express the scattering coefficients in terms of these solutions as 

t{x, z- k) = l/a+, Ri{x, z- k) = -P~/a+, (F.5) 

a+ = --^^T^ — ^^Uk + tp'ix, km - tp+(z, A;)le(i/2)[^(^)-^(^)l-^^ 

2[k^+p+{x,k)p-{x,k)]V V , ;JL f \ , j\ 

+[k + ip+{x, k)] [k - ip-{z, A;)]e(^/2)[^(^)-^(^)l-''" }, 

^ :^ri:^— 3T7^7T:^^^|[^ + ^r(^,fc)][fc + ^P+(^^ 
2[k^ -\-p+[x,k)p~{x,k)\ I 

-[k + ip+{x, k)][k + ip-{z, A;)]e(^/2)[^(^)-^(^)l-''" }, (F.6) 
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where 77^ = j^p'^{w,k)dw. (The symbols and /9~ are in accordance with the 
notation used in section III of [16].) If /(— 00) = +00, then V{—oo) — +00, and (F.5) 
gives the asymptotic behavior of r(x, z] k) SkS z ^ —00: 

]-V{z) + r]{z,k) 



t{x, z; k) = C{x, k) exp 



[1 + 0(1)], 



(F.7) 



k)dw=- j -j^^dw[l + 0(1)], (F.8) 

where zq is a constant. The generalized transmission coefficient f(x, W] k) is the 
transmission coefficient for the potential that has a jump at the right end point x (see 
[16]). So it is obvious that f has the same asymptotic form as (F.7). The expression for 
Ri is also the same as that for Ri. We have 



t{x, z; W; k) — Ci (x, W, k) exp 

ik 



--V{z)+v{z,k) 



[1 + 0(1)1, 



Ri{x,z;W;k) = -1 



[1 + 0(1)]. 



(F.9) 
(F.IO) 



From (F.8) we can see that r](z) — o(\z\) as 2; — > —00. If 1// e \ then r](—oo, k) is 
finite, and so we may let 77 = in (F.9) by including e^(^°°''^) in Ci. Since limfe_>.oP^ = 0, 

we have C{k = 0) = e^^^)/^^ Ci{k = 0) = 2e^/\ and r]{z, k = 0) = 0. 

The results for the case /(— 00) = —00 can be obtained in the same way. The 
asymptotic expressions for f and Ri are obtained by replacing V with —V in (F.9), and 
by changing the sign of the ffist term on the right-hand side of (F.IO). 

(2) /(-oo) = c. 

We consider the second-order differential equation 



'^'Mz) + f\z)i^{z) + ^'^'^ 



—'il){z) — ik'il){z) 
dz 



(F.ll) 



(This is the equation satisfied by the quantities a~^{x,z) and I3~^{x,z) defined in 
section III of [16], as functions of z with fixed x.) We set ip{z) = exp [J^ q{w,k) dw], 
and substitute into (F.ll). This yields the differential equation 



^q{z, k) - Q^q{z, k) + q\z, k) = f{z) - k' - Q^ik. 
dz f[z) f[z) 



(F.12) 



Suppose that k"^ ^ (?. (Since we are interested in the region of small A;, we need not 
consider the case k^ — c^.) Then equation (F.12) has a solution that tends to — Vc^ — A;^ 
as 2; — > —00. Let us define, in terms of this solution q{z^ k), 

p{z, k) = q{z, k) + Vc^ - k'^, s(z, k) = [ik - q(z, k)]/f(z), 



'-){x,z,k) = - I \p{w,k) + p{w, —k)]dw, 6{x,z,k) = 
2 Jz 

Q(x, z, k) = iV — k'^(x — z) + 6(x, z, k). 



2i 



[p(w, k) — p{w, —k)] dw, 
(F.13) 
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t{x, z; k) = 
Ri{x,z;k) 



s{x, k)s{x, —k) 



1 - s{x, k)s{z, -A;)e2^0(^.^.'=) 
-s{z, k) + s{x, A;)e*®(^'^''=) 
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(F.14) 
(F.15) 



1 - s{x, k)s{z, -k)e'^iQ{x,z,k) ■ 

We can show from (F.12) that p{z,k) — B{k)[f{z) — c][l + o(l)] as z — >■ — oo, where 
B{k) is some ^-independent quantity. Since p{z, k) is o(l) as ^ — >■ — oo, obviously both 
^{x,z,k) and 6{x. z.k) arc o{\z\) as z — >■ — oo. It is not difficult to show that, in fact, 
7(2:, —00, A;) is finite. Note also that lim2_^_oo s{z, k) = {ik + ^ — /c^) /c. As long as 
\/ c^ — k"^ has a nonzero real part, e*®^^'^'^-* vanishes in the limit z — )■ —00, and we obtain 

f(x, z; W; k) = C2{x, W, k) exp [Vc^ - k'^ z + A;)] [1 + o(l)] , (F.16) 

Ri{x, z] W; k) = ~ (ik + Vc^ - k"^^ + o(l), (F.17) 

where d{z,k) = 9{zQ,z,k) with a fixed constant 2^0. This 9{z,k) is o(|2;|) as 2; — >■ —00. 
The expressions for r and i?; are of the same forms, with W replaced by V{x). 

(3) /(-oo) = 0. 

Substituting ^{z) = exp [ikz + p{w, k) dw\ into (F.ll) yields the differential equation 



2ik 



p{z,k)+p\z,k) = f^z). 



(F.18) 



This equation has a solution that has the asymptotic form, as z — >■ —00, 
p{z,k) = ^f\z)[l + o{l)] 



p{z,k) 



1 



b + 2ik 



-f{z)[l + o{l)] 



if 


lim 

z—^—oc 




= 0, 


(F.19a) 


if 


lim 


/'(^) 




(F.19&) 


if 


lim 

z—^—00 


/'(^) 


= 00. 


(F.19c) 



P{z, k) = -^f{z)[V{z) - V{-oc)][l + 0(1)] 

(Recall that f'{z) is monotone for sufficiently large \z\, by our assumption. The right- 
hand side of (F.19b) is replaced by zf^{z)[l + o(l)] when b + 2ik — 0.) Let p stand for 
the solution specified by (F.19). Just hke (F.13), we define 



s{z,k) = -p{z,k)/f{z), 
9{x, z, k) : 

Then we have 
t(x, z] k) 



'y{x,z,k) = -J \p{w,k) + p{w,-k)] dw, 
\p{w, k) — p{w, —k)] dw. 



^ik{x—z)+i0{x,z,k)+'y{x,z,k) 



1 

2i 



1 — s{x, k)s{x, —k) 



Ri{x,z;k) = 



1 - s{x, k)s{z, -A;)e2»['=(^-^)+^(2^'^.fc)] 
-s{z, k) + s{x, k)e^m^'-^)+0i.^,-^,k)] 



1 - s(x, k)s{z, -k)e'^i[H=^-z)+e{x,z,k)] ■ 



(F.20) 

(F.21) 
(F.22) 
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Prom (F.19) it follows that s{z, ±k) = o(l) as z — )■ — oo. It is not difficult to see that 
7(x, — oo, k) is finite. So (F.21) and (F.21) yield the asymptotic forms as z — >■ — oo 

f(x, z; W; k) = Csix, W, k) exp [-ikz + ie(z, k)] [1 + o(l)], (F.23) 

Ri{x, z; W; k) = Ci{x, W, k) exp [-2ikz + 2ie{z, k)] + o(l), (F.24) 

where 9{z,k) = 9{zo,z,k) with a constant Zq. This 9{z,k) is o{\z\) as z ^ — oo. If 
G -Fq then 6{—oo,k) is finite, and we may let = in these expressions by 
including e*^(~°°''^) in C3. The expressions for r and i?; have the same forms as (F.23) 
and (F.24), with W replaced by V{x). 

Appendix G. The existence of Rr{x, —00; k) 

We defined the reflection coefficients for semi-infinite intervals as Rr{x, —oo;k) = 
limy^_oo Rr{x-,y\ k). (For ImA; = 0, the limit e J, of A; + is implied when necessary, 
as in (2.4).) Here we show that such a limit exists for ImA; > 0. (The existence of 
Ri{oo,y; k) can be shown in the same way.) We use the integral representation[18, 16] 

Rr(x, y; k) = r f{zy(x, z; k) dz (G.l) 

and the asymptotic form of t{x, z; k) as z —00 given in appendix F. For A; = 0, we 
have exactly Rr{x,y; k = 0) = tanh{[V(|/) — V{x)]/2} (see (E.l)). We can let y —00 
in this expression (see (3.7)). In the following, we assume that k ^ 0. 

First, we consider the case /(— 00) = +00. Substituting (F.7) into (G.l) gives 

Rr{x, y; k) = C^{x, k) f f{z)[l + h{z)]e-^^'^+^'>^'^^'^ dz, (G.2) 

Jy 

where h{z) = o(l) as 2; — >■ —00. Since f{z) = — {1/2) {d/ dz)V{z), 

r /(^)e-^(^)+2''(^) dz = -L-yiy)+My) + r ^/(^)g-v(.)+2,(.) ^ (G_3) 

Jy 2 Jy 

where C is independent of y. In this case, V{z) tends to +00 faster than \z\. From 
appendix F we know that r]{z) — o{\z\) and r]'{z) — o(l) as 2; — >■ —00. So the first term 
on the right-hand side of (G.3) vanishes, and the second term is convergent, as y — > —00. 
It is obvious that the part including h{z) in (G.2) is also convergent as y — )■ —00, since 
the integrand decays even faster. Therefore, the limit y — >■ —00 of Rr{x,y;k) exists, 
irrespective of whether Imk > or Imk — 0. The argument is the same for the case 
/(-oo) = -00. 

Next, let us consider the case /(— 00) = 0. We write (F.21) as 

j^^-i[kz-e{z)]+-f{z) 

^(^'"'^^= l-^(.)e-^-[^--^(^)l ^^-^^ 

explicitly writing only the dependence on z, with u{z) = s{z, —k)s{x, k)e'^^''^ and 
B = [1 — s{x, k)s{x, —k)]e'^^^ . If ImA; > 0, then r{x,z;k) falls off exponentially as 
z ^ —00, and it is easy to show that (G.l) has the limit y ^ —00. So, let us assume 
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that Im/c = 0. If Im/c = 0, the quantities 9{z) and 7(2;) (defined by (F.20)) are real- 
vaued functions of z, since p{z, —k) = [p{z, k)]*. By our assumption, f{z) and f'{z) are 
monotone for sufficiently large \z\. From (F.18) and (F.19) it follows that p{z, ±A;), and 
hence 7(2;), 6{z), 0'{z), Keu{z), lmu{z), are all monotone for sufficiently large \z\. As 
can be seen from appendix F, 



e{z)^o{\z\), e'{z)^o{l), j{z)^0{l), u{z) = o{l) asz^-00. (G.5) 

Wc divide the integral of (G.l) in two parts as — J^^ + with some yo, and 
substitute (G.4) into the second part. Then 



ryo 

L 



/yo \ 
f{z)e^^^'^ ^nM"(^)e-2("+i)*[*^^-^(^)l dz. 



(G.6) 



.n=0 



Here we have expressed as a power series in terms of u. This infinite series is 
convergent if —yo is chosen to be sufficiently large. Let us define 



yo 



Setting Z = z — [9{z)/k], we can write 



nB' 



rz=yo 

/ A(Z)e-2("+i)^'=^ dZ, A{Z) = 

J Z=V 



(G.7) 



(G.8) 



.-mz)/k] ■ 

Since f{z), 7 (-2), u{z), 9'{z) are all asymptotically monotone, there exists a number 
w such that A{Z) is monotone for z < w, and from (G.5) we see that A{Z) — >■ as 
z — > —00. Let us take yo < Then, using the second mean value theorem of integral 
calculus, we find 



iRe/J < 2 



B' 



;/(yo)e' 



27(2/0) 



k - e'iyo) 



|Im/„| < 2 



B' 



■f{yo)e 



27(2/0) 



k - e'iyo) 



(G.9) 



When y is finite, term-by-term integration is permissible on the right-hand side of (G.6). 
The infinite series J2'^=o^n is convergent and is equal to the left-hand side of (G.6). 
Taking the sum of (G.9) over n, we obtain 



Re 



yo 



f{z)T'^{x, z; k) dz 



< 2 



^2/(?/o)e 



27(2/0) 



k - e'iyo) 



1 



1 - Hyo)y 



(G.IO) 



and similarly for the imaginary part. (We are taking — yo to be sufficiently large so that 
|ii(yo)| < 1 and also k > e'iyo).) The right-hand side of (G.IO) is independent of y, and 
vanishes as yo —00. Therefore, the limit y —00 of (G.6) exists. Hence it is obvious 
that limj^^_oo Rrix, y; k) exists. 

The case /(— 00) = c 7^ can be studied in the same way, using (F.14) instead of 
(F.21). If Re Vc^ — k^ 7^ 0, then T(x,y;A;) falls off exponentially as y — >■ —00, and it is 
easy to show the existence of the hmit y — >■ —00 of (G.l). If k is real and k^ > c^, the 
right-hand side of (G.l) oscillates as y — > —00 and does not converge. In order to make 
the integral have a definite value, it is necessary to to assume that k has an infinitesimal 
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imaginary part ie. Then lim^j^o lim^^-^-oo of (G-l) exits, as can be seen from the fact that 
the following limit exists: 



(Apart from a constant factor, the integrand on the left-hand side of (G.ll) is the 
asymptotic form of f{z)T'^{x,z;k) as z — )■ — oo, where K = \Jk'^ — (9- and Z = 
z — [9{z)/K]. The remaining part of f{z)T'^{x,z;k) vanishes as 2; — )■ —00, and its 
integral converges more rapidly.) In this paper, however, we need not be concerned 
with the case k'^ > c^, since we are studying the small- A; expansion. The case k'^ = (? is 
not discussed here, since we do not need this case either. 
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